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NOTATION 


Amplitude of incoming wave 
Waterplane area . 

Added mass coefficients 

Damping coefficients 
Hydrostatic coefficients 
Exciting force 

Gravitational acceleration 
Imaginary unit 

Moments of inertia (roll, pitch and yaw) 
Longitudinal metacentric height 
Transverse metacentric height 
Wave number of incoming wave - 
Mass matrix 

Mass of the platform 


Moment of waterplane area 


Components of unit vector directed into the fluid 


Longitudinal location of center of flotation 
Longitudinal location of center of gravity 
Vertical location of center of gravity 
Water density 


Displacement of platform 


Angle of incoming wave in the xy-plane positive clockwise 


(present computer program notation) 


Angle of incoming wave in the xy-plane positive counterclockwise 


(ITTC notation) 
Frequency of wave encounter 


Wave frequency 
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Total velocity potential 
Potential of incoming wave 


Disturbance velocity potential due to the presence of the 
platforn 


Velocity potential due to motion of the platform with unit 
amplitude in each of six degrees of freedom 


Amplitude of motion in each of six degrees of freedom 
(surge, sway, heave, roll, pitch, and yaw) 


Diffraction potential 


Figures 5 through 16, the following notation is used: 
Wave period | 
Free-surface wave elevation 
Saige amplitude of motion 
Sway amplitude of motion 
Heave amplitude of motion 
Roll amplitude of motion 
Pitch amplitude of motion 
Yaw amplitude of motion 


Phase angle of motion (e=0: maximum positive motion occurs 
when wave crest is amidship.) 
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\, ABSTRACT 


The motions of a semi-submersible platform are computed for three 
different headings relative to the waves. This was done as part 
of a study established to compare various methods of seni- 
submersible platform motion prediction organized by the ITTC 
Ocean Engineering Committee. A strip method is applied in the 
numerical analysis. The results of this analysis are presented 
alone, that is, no comparisions are made with other theoretical 
methods or experimental data. This can be done when the ITTC 
Committee publishes the results of all participants of this 
study. The results of heave motion show a resonance when the 
wave period is about 3.35 seconds. _ : 


ADIMINSTRATIVE INFORMATION 


This study was performed under the Department Overhead Function. 
The Work Unit number is 1500-001. 


INTRODUCTION 


The 17th ITTC Ocean Engineering Committee, which met in the Autumn 1982 at 
Tokyo, decided to sponsor a comparative study on the motions of a semi- 
submersible. The primary purpose of this study is to compare existing com- 
putational methods and to validate their results. An opportunity was thus 
provided to further validate the computer program currently used by the Special 
Ships and Ocean Systems Dynamics Branch at the David Taylor Naval Ship Research 
and Development center. 

The platform configuration to be used for the study was selected by 
the ITTC Committee, and its description is given later in this report. 

The existing computer program was developed for semi-submersible platforms, 
with strip theory being applied as described in Reference le 


EQUATIONS OF MOTION 


The coordinate system, oxyz, is fixed at the midship section of the semi- 
submersible platform (see Figure 1). The oz-axis is directed vertically 
upward, and the oxy-plane is in the plane of the undisturbed free surface. 

The total velocity potential of the fluid in the presence of the platform 
ie expressed as 


O(x,y,z,t) = Rel(oy + dp)e7tut} (1) 


*References are listed on page 6. 
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where $7 is the potential of the incoming wave and is given as 


igA 
oy = -——exp[Kz + ikxcos8 - iKysin§] (2) 
Wo 


and $g is the disturbance velocity potential due to the presence of the 
platform. In Equation (2), A is the amplitude of the incoming wave, is its 
frequency, g the gravitational acceleration, 8 the angle of the incoming 

wave relative to the ox~axis (8 = 0 deg is following seas and 8 = 180 deg is 
head seas), and K = Wo2/g is the wave number. The disturbance potential 
consists of the following velocity potentials . 


op = 7 + LEyWy (3) 


where yy (j = 1, 2, ..-, 6) is the velocity potential arising from unit 
saplitule of platform in each of the six degrees of freedom, and 

bj (j = 1, 2, ..., 6) is the amplitude of motion in each of the six degrees of 
fre- iom. The diffraction potential is represented by 7. 


The potential v5 is determined as the solution of the Laplace equation 
with appropriate boundary conditions. To avoid the difficulty in solving a 
three-dimensional numerical problem, a strip theory or two-dimensional 
method is applied to obtain Vj- The fundamental assumptions of strip theory 
are discussed in Reference 2. 

Using strip theory, py is solved for j = 2,.3, 4, and 7 (sway, heave, 
roll motion and diffraction potential) at each section, and potentials for 
pitch and yaw are computed by simply multiplying the longitudinal distance of 
each section from the origin by the heave and sway potentials respectively. 
The surge potential, ~}, is assumed to be zero. 

In applying potential theory, it is assumed that the fluid is incompres- 
sible and inviscid, and that the flow is irrotational. The potential, Wy 
(j = 2, 3, 4, and 7) is determined as the solution of the following equation 
with the conditions specified. 


1. Laplace equation in the fluid domain 
a2y, a2y, 
+ = 0, for j = 2, 3, 4, 7 
dy2 az2 : 


2. the body boundary condition 


ayy 


= ~iwny, for j = 2, 3, 4 
én 
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oo for j=7 
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SY 
rs 
be 3. the linearized free-surface condition 
re av 
se Ky - = 0, on z=0 for j=2, 3, 4, 7 (6) 
oz 

ed 
ae The right hand side of Equation (5) is the normal velocity component at the 
ra body and the unit vector is directed into the fluid domain. All potential 
ap functions, a must satisfy the radiation condition for outgoing progressive 
Val waves at infinity, and become zero as z becomes negative infinity. 

. The solution of Equation (4) with the boundary conditions, Equations (5) 
-y and (6), is given in Reference 3. The hydrodynamic forces and moments are 
; obtained by integrating the pressure on the wetted surface of the body. The 
‘ ‘ details of the derivation of these forces and moments for heave and pitch 
8g motions are given in Reference 4. The forces and moments for other motions can 
SS be easily derived using a procedure similar to that described in Reference 4. 
ht If we let ay = &ye714t and equate the inertia force with the hydrodynamic 
iS force, the eauae ton of motion can be expressed as 
1. 
sy oe e 
a EL(mje + Ayn day + Byaoy + Cyeaj] = fy (7) 
Cy where ™jk is the mass matrix, Ai, the added mass, Byk the damping force, Chk 
4 the hydrostatic coefficients and f; the exciting force. The equations for Ajk, 
M3 Bjx, and fy are given in Reference 4. The mass matrix, mj,, is expressed as 
Zt 
re mj = M, forj=k=1, 2, 3 
. m4 = Ty 
a m55 = I5 
hs m6 = 16 (8) 
) = z= 
oe m24 = m2 = —Mzg 
*R 
La M15 = M5) = Mzg 
nnd 
os all other my, = 0 
‘ where M is the mass of the platform; 14, I5, and Ig roll, pitch and yaw moments 
ing of inertia, and zg is the vertical location of the center of the gravity with 
1 respect to free surface. The hydrostatic coefficent, Chk» is given as 
ie 
3 
io C33 = P&Awp 
= : C35 = C53 = ~08(Myp - XgAwp) 
ae C44 = AGMy (9) 
te C55 = AGM, + PgAwp(xXg - x,)2 
ao 
mt all other Cj, = 0 
Wy 
ae 
a) 
hs 
bya 3 
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where Ayp is the waterplane area, Myp its moment, Xy the longitudinal center 
of gravity, xq the longitudinal center of flotation, GMr the transverse 
metacentric height, Gi, the longitudinal metacentric height, and A the 
displacement of the platform. 


NUMERICAL METHOD AND RESULTS 


The model description is given in Figures 2, 3, and 4, and its principal 
dimensions are given in Table 1. In the numerical analysis, only 8 vertical 
columns and lower hulls are included in the computation. The transverse, 
longitudinal and diagonal columns are small compared to the vertical columns, 
and therefore, are excluded from the computation. 

The platform is divided into 29 sections in the longitudinal direction. 

At each section, the hydrodynamic forces (added mass, damping force and exciting 
force) are computed, and later these sectional forces are integrated along the 
length. There are 17 sections with lower hulls only and 12 sections with lower 
hulls and columns. The sections with lower hull only are treated as completely 
submerged sections and those with hull and column are treated as floating 
sections. The locations and types of sections are given in Table 2. 

The computations have been carried out for three different headings: 6 = 0 
(following), 8 = 315 (x = 45), and 6 = 270 (x = 90). ‘The heading angle of the 
incoming wave in the oxy-plane is taken positive counterclockwise by the 
ITTC Committee, while this angle is taken positive clockwise in the existing 
computer program. Therefore, the relationship between yx and 8 is: xy + B = 360. 
A step increment of wave period of 0.2 seconds was used for periods between 
0.5 and 10 seconds. Wherever there is resonance in the motion, smaller time 
time increments were used. Furthermore, it has been assumed that the water 
depth is infinite even though it was given as 3 m by the ITTC Committee. 

This assumption is necessary because the present computer program is applicable 
only for the deep water. The computer program can be easily extended to treat 
the case of finite water depth. 

The results are plotted in Figures 5 through 16. Since the experimental 
or numerical results produced by the other participants are not available at 
the present time, only the results computed by the author are shown in the figures. 
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Table 1 : Principal Dimensions 


sa 
‘ay Length of lower hull, L . 1.797 m 
Beam (total) 1.172 m 
iM Beam (one hull) 0.234 m 
Hull spacing between centerlines 0.938 m 
Draft 0.313 m 
; Center of gravity at the midship 
3 at the center line 
« 0.281 
Mt, Metacentric height 0.037 
0.046 
4 Radius of gyration 0.536 
a” 0.556 
st 0.634 
re: Displacement in fresh water 230.3 
Ban Water depth infinite 
y Wave direction 360, 315, 270 
x 0, 45, 90 
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os Table 2 : Locations of Sections 
J 
Z d 
«4 
2 ; 1 (AP) 
Neg 2 Lower aah 
3 ci 
4 ° 
5 Lower Hull + Column 
6 oe 
7 we 
mat 8 Lower Hull 
pea 9 Lid 
ms ; 
aes P Lower Hull 
3eu : iy 
aS . oe 
Ex Lower Hull 
sy o 
"y ‘ . Lower Hull 
mas , ry 
we in 4 oe 
Ba e Lower Hull 
3 e 
ar ; 
ian e Lower Hull 
a ey ry 
& i. o 
555 
nd ° Lower Hull 
r¢ pas we 
3 , : 
a 29 (FP) 
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FIGURE 8 - AMPLITUDE OF ROLL MOTION 
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FIGURE 16 ~ PHASE ANGLE OF YAW MOTION AT y=45 
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